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1. Introduction 

The transport in Josephson junctions with various different materials constituting the 
normal region has been a very active research field for decades now and continues to 
be one. Among the many interesting transport properties of Josephson junctions are 
the so-called subharmonic gap structure (SGS) and the excess current, both of which 
were accurately explained by the concept of multiple Andreev reflections (MAR). 
The MAR theory was first formulated for normal-superconducting (NS) interfaces 
by Blonder, Tinkham, and Klapwijk (BTK) [1, 2] and the BTK theory and its 
extensions (especially to the ferromagnetic or non-BCS superconducting contacts) 
are still actively used in fitting experiments [3, 4, 5] and in theoretical studies [6, 7]. 

The BTK theory was then extended to full SNS junctions by Octavio, Tinkham, 
Blonder and Klapwijk (OTBK) [8] and Flensberg, Bindslev Hansen and Octavio [9]. 
The OTKB approach does not keep track of the evolution of the quasiparticle phase 
between the interfaces and therefore assumes complete dephasing in the junction 
area. This assumption breaks down for sufficiently small systems such as, e.g., 
atomic wires and the fully coherent approach developed in mid 90's [10, 11, 12] is 
applicable instead. Nevertheless, the OTBK theory describes certain systems, such as 
microbridges, very well and keeps on being used in the literature both in experimental 
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[13, 14, 15, 16, 17, 18] as well as theoretical [19] studies. In particular, its extension 
to the experimentally relevant situation of asymmetric junctions was developed and 
applied in Refs. [13, 14, 15]. 

In this work we analytically study the excess current in the OTBK theory. 
Although the excess current has been derived analytically in more recent coherent 
theories [11] it has not been reported yet in an analytic form in the older incoherent 
OTBK approach. We fill in this gap and provide the analytical derivation of the 
excess current for incoherent generally asymmetric SNS junctions described by the 
OTBK theory. Our formula can be used for the experimental fitting but it also has 
implications for the understanding of the role of coherence within the junction as 
discussed in more detail in the concluding section. 

2. OTBK model 

The BTK theory describes the transport through a single normal-superconducting 
interface, which is assumed to consist of a ballistic superconductor in contact with an 
equally ballistic piece of normal metal. Scattering can thus only occur at the interface 
at x = 0, which is modelled by a repulsive delta-function potential HS(x), H > 0, 
with a dimensionless parameter Z = H/hvp (vf being the Fermi velocity) that 
represents the barrier strength [2]. Transport properties are found by matching the 
wave functions on either side of this barrier. The different allowed processes arc 
identified and labelled as follows: Andreev reflection A, normal reflection B, and 
transmission T. The corresponding probabilities A(E), B{E) and T(E) are expressed 
as functions of the quasiparticle energy E, the superconducting gap A, and the 
interface's barrier strength Z (cf. Table 1). The electrons at the normal side of 
the interface are separated into left- and right-moving populations, represented by 
the distribution functions f^(E) and /_>(£?), respectively. The current through the 
interface is then given by 

i= ik>rj E{jAE) - fAE)h (i) 

where Rq = (2N(0)e 2 VFA) 1 is the Sharvin resistance of the perfectly transparent 
interface (Z = 0) with A being the effective cross section of the contact and N(0) the 
(single spin) density of states at the Fermi energy Ep. Blonder, Tinkham and Klapwijk 
showed [2] that the distribution function for the left-moving electrons is given by 

.MS) = A(E) [1 - .M-£)] + B(E)f^(E) + T(E)f (E), (2) 

with fo(E) being the thermal Fermi distribution function fo(E) = 1/(1 + exp(/3(E — 
fj-))), assuming that the incoming electrons are in thermal equilibrium with their 
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Table 1. The reflection and transmission probabilities for an AfS-interface with 
the dimensionless barrier strength Z (after Ref. [2]; modified). 
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Figure 1. Schematic representation of an SNS junction, with arrows indicating 
the allowed processes at the interfaces: Andreev reflection A, normal reflection B 
and transmission T. 
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respective leads at the temperature l/ksP and the chemical potential \i. 

This description was extended in Ref. [8] to an SN-intcrface followed by an 
NS-intcrface, i.e. to a full Josephson junction (cf. Figure 1). We assume the same 
superconducting material on both sides, i.e. the same superconducting gap A, but 
different contacts and therefore differing barrier strenghts. We have interface 1, located 
at x = with barrier strength Z\, reflection and transmission probablities Ai(E), 
B 1 {E) and, T^E) and interface 2, at x = L with Z 2 , A 2 (E), B 2 (E) and T 2 (E). The 
distribution functions f+±(E, x), which are again to be taken in the normal region, are 
also functions of the longitudinal position within the junction, x. Now Eq. (2) can be 
applied to each of the two interfaces, which yields the following two equations [8] 

,M£,0) = A^E) [1 - f^(-E,0)] + B 1 (£7)/^(£7,0) +T 1 (E)f (E), (3) 

f^(E,L) = A 2 (E) [1 - f^(-E,L)} + B 2 (E)f^(E,L) + T 2 (E)f (E). (4) 

Note that we only combine distribution functions and not the quantum states 
(wavefunctions) at both interfaces. The relative phase of those states is therefore 
not considered, which is why the OTBK model only applies to incoherent junctions. 

Since all energies are measured with respect to the local chemical potential, right- 
moving quasiparticles with energy E at x — will arrive at x — L with energy E + eV, 
while left movers with energy E at x — L will have energy E — eV at x — 0. Thus the 
distribution functions at the interfaces relate to each other as 

f^(E,L) = f^(E-eV,0). (5) 

Eqs. (3)-(4) can be combined to eliminate, e.g., the left-moving part. Using Eq. (5) 
we can also shift all distribution functions from x — L to x — and hence omit the 
position argument in the following. The following equation can then be derived [13] 

/- (E)=A 1 (E){l-A 2 (-E + eV) [l - (E - 2eV)\ 

- B 2 (-E + eV)U(-E) - T 2 (-E + eV)f (-E + eV)} 

+ B^E^A^E + eV) [l - f^(-E - 2eV)] (6) 

+ B 2 (E + eV)f^(E) + T 2 (E + eV)f {E + eV)} 
+ T 1 (E)f (E) 
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that couples /_>(£) with /-(-£), - 2eV) and /_>(-£ - 2eV) and thus gives 

rise to an infinite system of linear equations for, say, f_,(E). 

3. Equal barriers 

The simplest far as the barriers are concerned, is the the case in which both 

interfaces are characterized by the same barrier strength Z\ = Z2 = Z. For this 
case an additional relation f^(E,0) = 1 — /<-(— E, L) was derived in Ref. [9] from 
Eqs. (3)-(4) and substituted into Eq. (5), which yields 

fc(E) = l-fe(-E-eV) (7) 

and greatly simplifies the problem. As before, the suppressed position arguments 
imply x = 0. Using this result we can reformulate Eq. (1) to depend on right-movers 
only 

1 = 4" r dE CM^) + f^~ E - eT/ ) - X ) • ( 8 ) 
eRo J -00 

Furthermore we can make use of Eq. (7) to eliminate the distribution functions for 
left-moving electrons in Eq. (3), which yields a significantly simpler equation than the 
fully general one from OTBK (6), namely 

yU(£) = A(E)f^(E - eV) + B{E) [1 - - eV)} + T(E)f (E). (9) 

The infinite system of linear equations generated by Eq. (9) was solved numerically 
in Ref. [9] to obtain subharmonic gap structure and excess current, but the latter can 
be obtained analytically [20], as we reproduce for convenience of the reader in the 
following. 

3.1. Normal current 

We shall first calculate the normal current to demonstrate the course of the derivation 
and to define some of the quantities used later on. We introduce the reflection and 
transmission probabilities in the normal case 

S » = IT^< Tn = 1 - Sn = TT^' (10) 

which are indeed the limits of B(E) and T(E) for vanishing A, as can be seen from 
Table l.f This allows us to rewrite Eq. (9) for the normal case as 

ft{E) - B n [1 - f^(-E-eV)]+T n f (E), (11) 

where (E) is the right-moving distribution function for the normal case. We rewrite 
the above Eq. (11) for the energy — E — eV 

ft(-E - eV) = B n [1 - f^(E)} + T n f (-E - eV), (12) 

insert this again in Eq. (11) and solve for f^(E), which yields 

/" ( E ) ^ T^fo{E) + ^kfo(E + eV) 

X In the normal state the Andreev reflection coefficient A(E) is identically zero. 
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where we have used that fo(—E — eV) = 1 — f (E + eV). We can write down the 
integrand from Eq. (8) with these normal-case distribution functions and simplify it 
to give 

-f (E+eV) 

ft(E) + ft(-E -eV)-l = ±±^(f (E) +'f (-E-eV)-i) 

+TT^( f °( E + eV ^ ) - l ) ( l4 ) 

-fo(-E) 

= {fo(E)-fo(E + eV))/(l + 2Z 2 ), 

which is easily integrated and yields the familiar result I n = V/R n , with the normal 
state resistance R n = (1 + 2Z 2 )R of the two-interface ballistic sandwich. § 

3.2. Excess current 

In the superconducting case we are interested in the excess current 7 0XC defined as 
J cxc = I — I n in the limit eV — > oo, which is what we will assume in the rest of this 
section. We define Af_>(E) = /-(£) - f^(E), insert this into Eq. (8) and substract 
the normal part thus arriving at the formula for the excess current 

i r 00 

/oxc = lfu> I dE ( Af{E) + Afi ~ E ~ eY) ^ (15) 

where we have dropped the arrows from the notation, as we are only dealing with 
right-movers in this section. To calculate Af(E) we define, similarly to the above, 
AB(E) = B(E) - B n and AT(E) = T{E) - T n , substitute all these definitions into 
Eq. (9), and solve for Af(E) to obtain 

Af(E) = A(E) [f n (E - eV) + Af(E - eV)] 

+ AB(E) [1 - U-E - eV) - Af(-E - eV)] (16) 

- B n Af(—E - eV) + AT(E)fo(E). 

Examining Eq. (16) and keeping in mind that A(E), AB(E) and AT(E) tend toward 
zero for \E\ 3> A, we see that there exists only a certain energy range e of the order 
of a few multiples of A where A(E), AB(E), and AT(E) can be considered non- 
zero such that f(E) will only differ significantly from f n (E) within e around E = 
and E = —eV. Since we only consider large bias those two energy regions are well 
separated. Therefore we can split Af(E) into one part which is only nonzero for 
\E\ < e and vanishes for all other energies and one part with the same properties for 
\E + eV\ < e. We introduce these parts by writing 

Af(E) - Af(E) + Af- eV (E). (17) 

This mathematical procedure is fully in line with the physical intuition that the only 
changes of the distribution functions induced by the superconductivity will occur 
within the few-A-multiplcs vicinity of the two Fermi energies of the leads. We can 
rewrite Eq. (16) for Af(E — eV) and we see that for \E\ < e most terms in the right 
hand side simply drop out as they include a vanishing multiplier. So we are left with 

Af~ eV (E - eV) = -B n Af(-E). (18) 

§ Note, that due to the ballistic nature of the junction this resistance is not just the sum of the two 
series resistances of the individual interfaces. 
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We insert Eq. (18) into Eq. (16), still assuming \E\ < e, to obtain 
Af(E) = A(E) [U(E - eV) - B n Af(-E)] 

+ AB(E) [1 - U-E - eV) + B n A.f(E)] (19) 

+ B 2 n Af(E) + AT(E)f (E). 

It turns out to be convenient to consider the combination Af(E) + Af(—E) in the 
following, i.e. to symmetrize the problem. Therefore, we rewrite Eq. (19) for Af(—E), 
sum the result with Eq. (19) and solve for Af(E) + Af(—E) obtaining 

1 + B n (A(E) - B(E))] (Af(E) + Af(-E)) = 

(A(E) - AB(E)) [f n (E - eV) + f n (-E - eV)} + 2AB(E) + AT(E). (20) 

Obviously, the above Eq. (20) again holds only for the range \E\ < e in which it was 
derived. The sum f n (E — eV) + f n (—E — eV) that turns up on the right-hand side 
of Eq. (20) can be calculated using Eq. (13) and the assumption of large bias, i.e. 
eV — > oo, so that we obtain 

>2 



f n (E - eV) + U-E - eV) = l±^'(f (E - eV) + f (-E - eV)j 

ITT^ P^), ^. (2D 
i 

From the condition of probability conservation A(E) +B(E) +T(E) = 1, we see that 
AT(E) = -A(E) - AB(E). We substitute this into Eq. (20) along with Eq. (21) and 
are left with 

(A(E)-AB(E))(1 + Z 2 ) 



Af(E) + Af(-E) 



(22) 



l + B n (A(E)- B(£))j(l + 2Z 2 ) 
We now take the integrand from Eq. (15) and expand it by inserting Eq. (17) to 

get 

Af(E)+Af(-E-eV) = Af(E) + Af- eV (E) + Af(-E-eV) + Af- eV (-E-eV).(23) 

The first and the last terms are nonzero around E = and we can use the 
straightforward modification of Eq. (18) for the simplification Af(E) + Af~ eV (—E — 
eV) = (1 — B n )Af(E). Analogously, the two middle terms in Eq. (23) are nonzero 
around E — —eV and since they appear under the integral extending over the entire 
energy range and have strongly localized support their energy arguments can be shifted 
so that they are localized around E = as well.|| This eventually leads to the relations 
Af- eV (E) + Af(-E-eV) Af- eV (E-eV) + Af(-E) = (l-B n )Af(-E). Putting 
all the pieces together leaves us with 

/oxc = ^ j°° dE (1 - B n ) (A/(£) + Af(-E)) 

= 1 r rlF ME)-AB{E) 

eR (l + 2Z*)J_ oo [l + B n (A(E)-B(E))]' 1 ; 

where we used Eq. (22) to produce the final integrand. The analytical integration must 
be performed for \E\ < A and \E\ > A, separately, because A(E) and B(E) take on 

|| The strongly localized support of the involved terms is essential for the possibility of the variable 
shift and its lack can lead to seemingly paradoxical results when done formally, e.g., in Eq. (8). 
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Figure 2. The excess current of a symmetric, fully incoherent Josephson junction 
in the OTBK model as a function of the barrier strength Z. 



different functional forms in these intervals. It can be evaluated using trigonometric 
or hyperbolic substitution for the subgap or overgap energies, respectively. Thus we 
find the excess current in the symmetric case is given by 

e/c ^ fln = 2(1 + 2Z 2 )tanh- 1 (2Zy/(l + Z 2 )/(l + 6Z 2 + 4Z 4 )) 

x (V(1 + Z 2 )(1 + 6Z 2 + 4Z 4 )) ~* - ~. (25) 

The first (and longer) term on the right hand side of Eq. (25) results from the subgap 
integral, the —4/3 term is the contribution of the overgap part. This analytic result is 
plotted in Figure 2 and the comparison with the earlier numerical result by Flensberg, 
Bindslcv Hansen and Octavio plotted in Figure 6 of Ref. [9] shows a nice agreement. 

4. Differing barrier strengths 

To obtain a similar expression for the excess current in the case of asymmetric barrier 
strenghts we need to restart from Eq. (6), since Eq. (7) and the ensuing simplifications, 
in particular Eq. (9), cannot be used. The course of the derivation, however, is very 
similar to the above. Again we start by calculating the normal current. 

4-1. Normal current 

Corresponding to Eq. (1), the normal current is given by 

In = 4" r dE (/ " {E) - f ^ E)) ' (26) 
e-Ko J-oo 

where f^(E) is defined as above and f^_(E) is its left-moving counterpart. The 
reflection and transmission probabilities in the normal case, B n i and T a ^, are defined 
as above with the additional index i g {1; 2}, which indicates the interface in question. 
From Eq. (6) we now find for the right-movers in the normal case 

ft(E) = B nA [B n , 2 f^(E) + T na f (E + eV)] + T nA f (E). (27) 
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Since we also need to consider left-movers in this section, we use Eqs. (3)-(5) to derive 
the left-moving counterpart to Eq. (6), which is not shown for reasons of length, and 
finally the equivalent of the above Eq. (27) for left movers, which reads 

ft(E) = B n ,2 [B nA ft(E) + T nA f (E)} + T n , 2 / (£ + eV). (28) 

We take the integrand from Eq. (26) and use Eqs. (27)-(28) to rewrite it as follows 

fZ(E) fl(E) = Tn f n ^ (f (E) f a (E + eV)) 

= 1 + z l + Z 2 (fo( E )-ME + eV)), (29) 

This is easy to integrate and yields I n = V/(p n Ro), with p n = 1 + Z\ + Z\. Note that 
p n Ro simply becomes R n for Z\ = Z 2 = Z, so we find the normal current from above 
for equal barriers again. 



4-2. Excess current 

For the calculation of the excess current we assume large bias once again and introduce 
AM(E) = M(E) - M(E), AB t (E) = B t (E) - B n<i and AT,(£) - Ti(E)-T ati , just 
like above in the case of symmetric barriers. Using these relations we can expand 
Eq. (6) and subtract Eq. (27) to obtain 

AM {E) = A 1 {E){l-A 2 (-E + eV)[l- M (E - 2eV)] 
-B 2 (-E + eV)M(-E) 
-T 2 (-E + eV)M-E + eV)} 

+ B nA {A 2 {E + eV)[l - M{-E - 2eV)] + B a , 2 AM(E) (30) 
+ AB 2 (E + eV)M(E) + AT 2 (E + eV)f (E + eV)} 

+ ABi(E)^A 2 (E + eV) [l - M(~E - 2eV)] 

+ B 2 (E + eV)M(E) + T 2 (E + eV)f (E + eV)} 

+ AT 1 (E)ME). 

By the same logic as before we see that Af^(E) is only nonzero for \E\ < e or 
\E + eV\ < e. Therefore we split Af^(E) into two parts, just like we did above and 
with the same properties 

AM(E) = AM(E) + AM eV (E). (31) 

For the remainder of the section we assume small energies (\E\ < e), in which case 
Eq. (30) can be reduced and solved for A /_>(£?) to yield 

[1 - Bi(£?)B n>2 ] AM(E) = A 1 (E)[l - B na (M{-E) + AM(-E)) 

~M 2 M-E + eV)] (32) 

+ AB 1 (E) [B nt2 M(E) + T n , 2 / (E + eV)} 

+ AT 1 (E)ME) 
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We rewrite Eq. (32) for Af_>(—E), sum the result with Eq. (32) and solve for 
A/_> (E) + Af^(-E), which gives 

[1 + B^A^E) - B^E))} [A/_> (E) + Af^(-E)] 

= A 1 {E)[2-B n>2 (ft (E) +ft(~E)) 

-T», 2 (/o(-B + eV) + f Q (-E + eV))] 
+AB 1 (E)[B n , 2 {fZ(E) + ft{-E)) (33) 

+TnAfo(E + eV) + fo(-E + eV))} 
+AT 1 (E) (fo(E) + fo(-E)) 

" v ' 

1 

The sum fo(E + eV) + fo(—E + eV) in the above becomes zero for large bias, which 
means that the terms explicitly involving T n , 2 drop out of Eq. (33). Furthermore, 
using Eq. (27) we can write 



ft(E) + fM-E) = [f? n ,iT n , 2 (h(E + eV) + f Q (-E + eV)) 

+T„,i (fo(E) + /o(-25))]/(l - B n sB na ) = ^ , (34) 

v ' 1 — -£>n,l-Dn,2 

1 

further simplifying Eq. (33), which can now be written as 

Pn 1 + 5 n , 2 (^l(£) - 

where pi — 1/T n ^ = 1 + Z"f is the dimensionless resistance of the single i-th interface 
in the normal state. In a similar way and using the same assumptions, i.e. eV — > oo 
and \E\ < e, we can show that 

Af^(E - eV) + A/--(-£ - eV) = - P j- - ^ a( f ~ ^ a( f fr, tl . (36) 

Pn 1 + Bn4(A 2 (S) - B 2 (-E)J 

We still need to get the left-moving equivalents of Eqs. (35), (36), so first we 
derive the counterpart to Eq. (30) for the left-movers, which is not shown, because 
the derivation follows the earlier pattern and does not deliver new insights. Just like 
above we can split Af<-(E) up into 

A/+_ (E) = A./L (E) + Af- eV (E) . (37) 

As for the right-movers and in much the same way we can show that for \E\ < e and 

eV — > oo 

A/-(*) + AU-E) - ^ i-^MaA (38) 

Pn 1 + Bn, 2 (-4l(^) - Bi{E)) 

as well as 

- .V, + A/Z-(- E - «V) - -g 1 + ff A - ( ^ (E)) - (») 



The excess current is now given by 

1 r°° 
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(35) (36)' 

A. A. 

d£7( A/^(S) + A/_(-.E) + Af- eV (E) + Af- eV (-E) (40) 



- [A/L(i5) + A.M-£)] - [A/Z eV (25) + Af- eV (-E)}). 

V v ' S v ^ 

(38) (39)' 

The braces and brackets in Eq. (40) indicate which terms in the expression correspond 
to which one of the above equations. The primed brackets are shifted in energy, which 
does not matter to the final result, since the integral extends over the entire energy 
range and the integrands have strongly localized support. Finally we can express the 
excess current as 



AAE)-ABAE) , x , A 



j_ r Am mb>- 

eRoJ-oo [Pnl + S n , 2 (Ai 

The integral in Eq. (41) can be solved and the result for Z\ > Z 2 is given by 
£lc ^ Ro = 2 Pn tanh- 1 (2Z^ Pl /(2Zl + (l + 2Z^ 

x (z iy / Pl (2Zl + (l + 2Z*)*) S ) 

l y/(Z* Zi)/p D tanh- 1 y / (Z 1 2 - Zl)/p n 



tan 1 



(42) 



x (1 + 2Z 1 2 )(1 + 2Z|) 1 
2JpZ{Zl - zl)i 

For Z\ < Z 2 the excess current is obtained by exchanging Z\ and Z 2 in Eq. (42) and 
the result is thus symmetric with respect to the interchange of the two interfaces. 
Again, the first term (the first line) on the right hand side of Eq. (42) corresponds 
to the sub-gap integral and it is easy to see how for Z\ = Z 2 = Z it becomes 
the corresponding term in Eq. (25). The remaining two terms (the second line), 
which result from the over-gap integral, converge towards —4/3 for Z\ — > Z 2 , as 
we now show by the Taylor expansion of tan _1 (z) = z — + — ^z 7 + . . . and 
tanh _1 (z) = z + iz 3 + iz 5 + iz 7 + . . . resulting in tan" 1 (z)-tanh _1 (z) = -|z 3 + (9(z 7 ) 

3 

so that the square bracket in the second line of Eq. (42) tends to — |((Z 2 — Z|)//0 n ) 2 
for Zi -» Z 2 . Therefore the last two lines of Eq. (42) reduce to _Ii± 2 £^±H£ll _ i j 
which simply becomes —4/3 for Z\ = Z 2 and we thus recover Eq. (25) in the symmetric 
case. The full result (42) is plotted in Figure 3 as a function of the two barrier 
strengths Zi 2 - The negative excess current predicted for large enough normal-state 
resistance persists for asymmetric junctions Z\ ^ Z 2 with arbitrarily large asymmetry 
although its magnitude decreases (also note the prefactor p n customarily multiplying 
the plotted excess current) and, thus, its experimental observation may be impeded 
by the asymmetry of real junctions. 



5. Conclusions and outlook 



In this work we have analytically calculated the excess current within the OTBK 
theory describing fully incoherent SNS junctions. We have confirmed previous 
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Figure 3. The excess current of an asymmetric, fully incoherent Josephson 
junction as a function of the two barrier strengths Z\ and Z2. The isoline 7 eX c = 
is shown in red. 

numerical findings [9] of negative excess current for large enough normal-state 
resistance in junctions with symmetric barriers. Furthermore, we extended those 
calculations also to the case of asymmetric barriers with qualitatively similar results, 
i.e. occurrence of negative excess current regardless of the asymmetry. Our formula 
(41) can be used also in the most general case of different superconducting leads, for 
an experiment see, e.g. Ref. [15], where Ai ^ A2. The presence of two gap values 
prohibits further analytical treatment, however, Eq. (41) still holds and the integral 
can be easily evaluated numerically. 

The numerical findings of Ref. [9] were challenged in Ref. [11] (p. 7372, paragraph 
below Eq. (30)) and the negativity of the excess current was interpreted as possibly 
stemming from a lack of convergence of the numerical study, i.e. from not reaching 
the true eV — ► 00 limit. Our study clearly demonstrates that this objection cannot 
hold since we explicitly work in the required limit, thus avoiding any finite-^ issues. 
We, however, do not question the presence of non-trivial issues in the experimental 
determination of the excess current related to the finite voltage and possible heating 
effects nicely reviewed and discussed in Ref. [12]. Apparently, observations of negative 
excess current (so called deficit current) have been reported in experiments [13, 14]. 

Nevertheless, we analytically prove the correctness of the old numerical results 
[9] predicting negative excess current within the OTBK theory. The discrepancy with 
the results of Ref. [11] then must stem from the difference of the two considered 
models, more specifically, the role of internal coherence of the junction. While the 
OTBK theory only considers matching of the distribution functions between the 
two interfaces corresponding to fully incoherent junctions, the Hamiltonian theory of 
Ref. [11] matches the wavefunctions throughout the whole junction thus fully retaining 
the coherence within the junction. The high- voltage properties of the two models differ 
even qualitatively, one predicting negative excess current for small transparencies, the 
other one not. Another qualitative difference between OTBK and the fully coherent 
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theory is in their dependence on the junction asymmetry: While the fully coherent 
results in the limit of strong coherent coupling to the leads (r 12 S> A, relevant for 
many experiments, e.g. [21, 22, 23]) only depend on the asymmetry through the total 
junction resistance [11], it is not so in the OTBK case as we can immediately see from 
our result for the excess current (Eq. (42)) which is not a function of Z\ + Z\ only. 

This finding shows that the coherence within the junction plays a crucial role 
for the superconducting transport even at finite voltage bias and therefore the level 
of decoherence/dephasing within a junction should be carefully considered when 
describing a particular experiment. This effect, i.e. nonzero dephasing within the 
junction, may be responsible for the experimentally observed discrepancies between 
the experiments [21, 22, 23] and theoretical predictions [11] systematically reported 
recently in Josephson junctions made of carbon nanotubes. While those discrepancies 
are currently interpreted as the superconducting gap renormalization this picture does 
not seem to be fully consistent with the positions of the subharmonic gap structure 
features, which appear at the positions determined by the un-renormalized gap value. 
The dephasing picture could capture the relevant physical mechanism instead although 
this remains an open issue in the currently booming field of superconducting transport 
in carbon-allotropes-based Josephson junctions. 

Apart from the obvious usage of our newly derived analytical formulae for the 
OTBK excess current to fit experiments for relatively large and thus fully incoherent 
junctions, they can also be used as a limit benchmark of future partially-coherent 
theories, possibly relevant for current nanoscale experiments. These experiments 
as well as future devices built from novel low dimensional materials with peculiar 
electronic structures, such as graphene nanoribbons, could realistically be described by 
existing dephasing approaches for atomistic models [24, 25] coupled to non-equilibrium 
transport. The development of such a partially-coherent theory and the analytical 
evaluation of its excess current interpolating between the two limits is our next step. 
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